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Abstract 


The  potential  for  using  a  digital  automatic  gain  control  (DAGC)  is  examined.  The  DAGC  operates  by 
digitally  storing  a  clutter  map  of  the  average  of  past  returns.  This  average  clutter  map  is  digitally  filtered  to 
remove  small  spatial  variations  and  the  resulting  values  are  used  to  adjust  the  front  end  receiver  gain.  In 
contrast,  a  traditional  analog  automatic  gain  control  (AAGC)  uses  real  time  filtered  feedback  to  control  the 
adjustable  front  end  gain.  The  physical  realizability  constraints  of  an  AAGC  typically  result  in  a  wider 
shadow  zone  for  regions  where  the  scan  is  moving  from  a  land  to  sea  clutter  environment.  It  is  shown  that 
the  DAGC  filtering  operation  on  the  clutter  map  of  past  returns  allows  the  spatial  extent  of  the  filter  to  be 
distributed  evenly  across  the  land-sea  transition  regions  thereby  reducing  the  shadow  zone.  The  AAGC  is 
shown  to  display  slightly  superior  noise  performance  over  the  DAGC  approach.  Simulations  of  the 
behaviour  of  both  systems  under  simplified  configurations  are  also  shown  to  qualitatively  illustrate  and 
clarify  the  tradeoffs  associated  with  either  approach. 

Resume  analytique 


On  etudie  la  possibility  d’utiliser  une  commande  numerique  de  gain  automatique  (CNGA).  La  CNGA 
fonctionne  en  stockant  une  carte  de  la  moyenne  des  precedents  retours  de  fouillis.  Cette  carte  de  fouillis  est 
filtree  numeriquement  afin  d’en  eliminer  les  petites  variations  spatiales.  On  se  sert  alors  des  resultats 
obtenus  pour  regler  le  gain  recepteur  avant.  Par  contraste,  une  commande  analogique  de  gain  automatique 
(CAGA)  classique  utilise  une  retroaction  filtree  a  temps  reel  pour  regler  le  gain  avant.  Les  contraintes  de 
realisabilite  physique  d’une  CAGA  produisent  generalement  une  zone  d’ombre  plus  grande  pour  les 
regions  ou  le  balayage  passe  d’un  environnement  de  fouillis  de  sol  a  un  environnement  de  fouillis  d’ocean. 
On  remarque  que  V operation  de  filtrage  de  la  CAGA  effectude  sur  la  carte  des  precedents  retours  de  fouillis 
permet  de  distribuer  uniformement  l’etendue  spatiale  du  filtre  d’un  cote  a  l’autre  des  regions  de  passage 
terre-mer,  reduisant  ainsi  la  zone  d’ombre.  II  se  revele  que  la  CAGA  presente  un  niveau  de  bruit  legerement 
superieur  a  celui  de  la  CNGA.  Des  simulations  du  comportement  des  deux  systemes,  accomplies  selon  des 
configurations  simplifies,  sont  egalement  montrees  pour  illustrer  et  clarifier  qualitativement  les  avantages 
et  les  inconvenients  lies  a  chaque  option. 
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Executive  Summary 


The  development  of  the  digital  scan  converter  (DSC)  and  its  inherent  digital  processing  capabilities  opens 
the  door  to  a  variety  of  different  schemes  for  digitally  manipulating  scanned  data  so  as  to  improve  the 
overall  detection  capability  of  the  AN/APS-506  surveillance  RADAR.  This  report  examines  one  possible 
technique,  namely  the  implementation  of  a  digital  automatic  gain  control,  DAGC,  as  part  of  the  RADAR 
system. 

The  DAGC  controls  the  gain  by  building  a  clutter  map  of  the  average  of  the  return  strengths  over  time.  The 
clutter  map  is  filtered  to  remove  small  spatial  variations  and  the  resulting  map  of  large  scale  variations  is 
used  to  predict  the  required  attenuation  for  future  incoming  returns.  Traditionally,  RADAR  systems  have 
employed  an  analog  automatic  gain  control  system  (AAGC),  which  is  essentially  a  ‘hard  wired’  system. 
Because  of  this,  it  is  difficult  to  adjust  an  AAGC  system  during  actual  operation  so  as  to  optimize  its 
performance  under  different  operating  modes.  This  often  results  in  overall  performance  degradation.  The 
digital  approach  offers  an  advantage  in  this  regard  since  its  processing  is  done  digitally  and  operating 
parameters  can  be  easily  adjusted  during  operation  to  allow  optimal  performance. 

This  report  focuses  on  the  development  of  the  operating  theory  for  both  the  AAGC  and  the  DAGC.  The 
stability,  frequency  response  and  noise  characteristics  of  both  systems  are  examined  and  compared.  In 
addition,  simple  simulations  are  performed  to  clarify  the  tradeoffs  associated  with  each  system.  It  is 
shown  that  both  systems  are  stable  over  the  full  range  of  realistic  operating  conditions.  The  DAGC  system 
is  shown  to  have  the  most  flexible  frequency  response  since  the  filtering  step  is  performed  digitally  and  the 
filter  shape  can  be  chosen  at  will.  This  flexibility  allows  the  DAGC  to  be  interactively  adjusted  to  provide 
a  finer  gain  control  resolution  at  the  land-sea  interface  thereby  reducing  the  spatial  extent  of  the  shadow 
region  that  the  ‘memory’  of  the  strong  land  clutter  returns  can  cause  over  the  sea.  The  practical  inability  of 
adjusting  the  RC  values  of  the  range  gated  filters  in  the  AAGC  usually  results  in  a  relatively  large  shadow 
region  along  the  coast.  While  the  AAGC  system  displays  a  small  superiority  in  noise  performance,  the 
difference  is  small  enough  that  it  is  not  likely  to  be  a  significant  factor  influencing  the  choice  of  system 
approach. 
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Sommaire 


La  mise  au  point  du  convertisseur  de  balayage  numerique  (CBN)  et  de  ses  inherentes  capacites  de 
traitement  numerique  ouvre  la  voie  a  tout  un  6ventail  de  moyens  qui  permettent  de  traiter  numeriquement 
des  donn6es  de  balayage  pour  ameliorer  la  capacite  de  detection  generate  du  radar  de  surveillance 
AN/APS-506.  Le  present  rapport  examine  une  technique  possible,  soit  la  mise  en  oeuvre  de  la  commande 
numerique  de  gain  automatique,  CNGA,  integree  au  systeme  radar. 

La  CNGA  regie  le  gain  en  eiaborant  une  carte  de  fouillis  repr&entant  la  moyenne  des  forces  des  retours 
enregistrees  au  fil  du  temps.  La  carte  est  filtr6e  pour  eliminer  les  petites  variations  spatiales  et  le  r^sultat, 
une  carte  aux  grandes  variations  d’echelle,  sert  k  prevoir  l’att6nuation  requise  pour  les  prochains  retours. 
Jusqu’ici,  les  systemes  radars  ont  utilise  un  systeme  de  commande  analogique  de  gain  automatique 
(CAGA),  un  systeme  essentiellement  «  cable  ».  De  ce  fait,  il  est  difficile  de  regler  le  systeme  CAGA 
pendant  1’operation  proprement  dite  pour  maximiser  sa  performance  dans  des  modes  de  fonctionnement 
diffcrents.  Cela  entraine  souvent  une  degradation  de  performance  generale.  Le  proced6  numerique  offre  un 
avantage  k  cet  egard,  etant  donne  que  son  traitement  se  fait  numeriquement  et  qu’on  peut  rdgler  facilement 
les  parametres  d’ exploitation  pendant  qu’il  fonctionne  pour  obtenir  une  performance  optimale. 

Le  present  rapport  se  focalise  sur  le  d£veloppement  de  la  theorie  d ’exploitation  visant  les  deux  systemes 
CAGA  et  CNGA.  La  stabilite,  la  r£ponse  en  frequence  et  les  caracteristiques  de  bruit  des  deux  systemes 
sont  examinees  et  compares .  De  plus,  de  simples  simulations  sont  effectu6es  pour  prdciser  les  a  vantages  et 
les  inconvenients  associ^s  a  chaque  systeme.  On  note  que  les  deux  systemes  s’av&rent  stables  dans  toutes 
les  conditions  pratiques  d’exploitation.  II  se  revele  que  la  r£ponse  en  frequence  du  systeme  CNGA  est  la 
plus  souple  puisque  l’etape  de  filtrage  se  fait  numeriquement  et  qu’on  peut  choisir  la  forme  du  filtre  k 
volont6.  Cette  souplesse  permet  de  regler  la  CNGA  de  maniere  interactive  pour  obtenir  une  plus  haute 
resolution  de  la  commande  de  gain  a  1’interface  terre-mer,  ce  qui  r£duit  l’etendue  spatiale  de  la  region 
d’ombre  que  la  «  memoire  »  des  puissants  retours  de  fouillis  de  terre  peut  provoquer  sur  la  zone  maritime. 
Puisqu’il  est  impossible  de  pratiquement  rdgler  les  valeurs  de  la  resistance  de  capacite  des  filtres  k 
ddclenchement  de  la  CAGA,  on  obtient  generalement  une  grande  rdgion  d’ombre  le  long  de  la  cote.  Si  le 
systeme  CAGA  presente  une  efficacitd  legerement  sup6rieure  quant  au  bruit,  la  difference  n’est  pas  assez 
importante  pour  etre  un  facteur  determinant  notable  dans  le  choix  de  systeme. 


McDonald,  Michael  K.,  Digital  Automatic  Gain  Control  (DADC)  for  Surveillance  Radar  Applications: 
Theory  and  Simulation.  ,  Centre  de  recherches  pour  la  defense  Ottawa,  DREO  TR  2000-115,  december 
2000.  (en  anglais) 


iv 


Table  of  Contents 


Abstract . i 

Resume  analytique . i 

Executive  Summary . iii 

Sommaire . iv 

Table  of  Contents . v 

List  of  Figures . vi 

List  of  Tables . vii 

1.0  Introduction . 1 

1.1  Receiver  Specifications . 3 

2.0  Automatic  Gain  Control  Theory . 4 

2.1  Analog  Automatic  Gain  Control  (AAGC) . 4 

2.2  Digital  Automatic  Gain  Control  (DAGC) . 9 

2.2.1  Transfer  Function  and  Stability  of  DAGC  without  Spatial  Filtering . 10 

2.2.2  Noise  Performance  of  DAGC . 13 

2.2.3  Frequency  Responses  of  DAGC . 16 

2.2.4  Transfer  Function  and  Stability  of  DAGC  with  Spatial  Filtering . 17 

3.0  Simulations . 22 

3.1  Configuration  1  Simulations  -  RADAR  Sweeping  from  Sea  to  Land  to  Sea . 25 

3.2  Configuration  2  Simulations  -  RADAR  Sweeping  Parallel  to  Land . 32 

4.0  Summary  and  Future  Work . 34 

References . 35 

List  of  Acronyms . 36 


v 


List  of  Figures 


Figure  1:  Typical  analog  AGC  loop  implementation . 1 

Figure  2:  Proposed  configuration  of  AGC  loop . 2 

Figure  3:  AGC  configuration  assumed  in  report . 5 

Figure  4:  Gain  control  characteristic . 7 

Figure  5:  Steady  state  static  regulation  curves  for  u  1=10,150  and  500 . 7 

Figure  6:  Sampling  action  of  range  gated  filter  banks.  Dotted  line  shows  impulse  response  of  RC  filter  in 
time  and  dashed  lines  show  how  impulse  response  is  distributed  in  time  by  range  gating  action.  Solid 
line  is  a  sample  sinusoid  input  signal.  Current  output  of  RGF  (TIME  =  0)  would  equal  integrated  sum 

of  dashed  filter  response  multiplied  by  sample  signal  from  TIME=0  through  to  (-  Infinity) . 8 

Figure  7:  Bode  plot  of  range  gated  filter  bank  and  AAGC  loop  with  feedback  arm  gain,  equal  to  150. 
Curves  are  displayed  for  average  AAGC  input  signal  power  levels  after  input  amplification  and  are  as 

follows:  receiver  noise  =  .025  W,  sea  clutter  =  0.65  W  and  land  clutter  =6.5  W . 9 

Figure  8:  Flow  chart  of  simplified  DAGC  loop  without  spatial  filtering.  N  equals  the  number  of  temporal 

samples  averaged  together  and  REF  is  the  desired  reference  level . 10 

Figure  9:  Noise  enhancement  factor  versus  accuracy  with  which  correction  can  be  applied.  A=1  represents 
perfect  application  of  desired  gain  correction  while  A<1  represents  an  underestimate  and  A>1  an 
overestimate  of  the  required  correction.  Curves  are  shown  for  cases  in  which  1,  3  or  10  past  scans  are 

averaged  to  form  the  gain  correction  estimate . 14 

Figure  10:  Step  response  of  simple  DAGC  loop  without  spatial  filtering  for  N=l,  3  and  10  scan  averaging. 


. 15 

Figure  11:  Bode  plot  of  DAGC  response  for  N=l,  3  and  10  scan  averaging,  x-axis  is  actually  w*T,  to  get 

frequency  divide  the  x-axis  values  by  the  sampling  period,  T . 16 

Figure  12:  Flow  chart  of  simplified  DAGC  loop  with  spatial  filtering.  See  text  for  description  of  variables 

and  subscripts.  See  figure  8  for  loop  definitions . 18 

Figure  13:  Test  configuration  1  to  simulate  the  effect  of  a  radar  beam  sweeping  from  sea  clutter  onto  land 

and  back  onto  sea  clutter.  Azimuthal  bin  width  is  79m  and  range  bin  width  is  75m . 23 

Figure  14:  Test  configuration  2  to  simulate  the  effect  of  radar  beam  sweeping  parallel  to  shore  with  AAGC 
sampling  window  of  3.8  ps  straddling  land-sea.  Azimuthal  bin  width  is  79m  and  range  bin  width  is 

75m . 24 

Figure  15:  Flow  chart  for  algorithm  corresponding  to  AAGC  processing  module . 26 

Figure  16:  Flow  chart  for  algorithm  corresponding  to  DAGC  processing  module . 27 

Figure  17:  Output  from  AAGC  for  configuration  1  deterministic  input  signal.  Arrows  identify  targets . 28 

Figure  18:  Output  from  AAGC  for  configuration  1  with  non-deterministic  input  signal  added.  Arrows 

identify  targets . 28 

Figure  19:  Plot  along  azimuthal  vector  containing  targets  of  DAGC  output  for  configuration  1  with 

deterministic  input  data.  Arrows  indicate  targets . 29 

Figure  20:  Output  from  AAGC  for  configuration  1  deterministic  input  signal  with  capacitance  reduced  by 

two  orders  of  magnitude  in  RC  filter.  Arrows  indicate  targets . 30 

Figure  21:  Plot  along  azimuthal  vector  containing  targets  of  DAGC  output  for  configuration  1  input  data 

with  non-deterministic  portion  of  clutter  signal  added.  Arrows  indicate  targets . 31 

Figure  22:  Plot  along  range  vector  containing  targets  of  AAGC  output  for  configuration  2  input  data  with 

non-deterministic  portion  of  clutter  signal  added.  Arrows  indicate  targets . 33 

Figure  23:  Plot  along  range  vector  containing  targets  of  DAGC  output  for  configuration  2  input  data  with 
non-deterministic  portion  of  clutter  signal  added.  Arrows  indicate  targets . 33 


vi 


List  of  Tables 

Table  1.1:  Receiver  Specifications . 3 

Table  2.1:  Regulation  Characteristics . 6 


1.0  Introduction 


,ni,erera  ran8e  deptndence  °f  RADAR  «*pW  with  .he  diverse  elude,  environmenls 

encountered,  means  .ha,  RADAR  signal  return  strengths  will  commonly  span  seven,!  otders  of  magnitude 

,nce  the  dynamic  range  of  the  receivers  is  typically  much  smaller  i,  is  necessary  ,o  provide  some  form  of 
gam  connol  ,o  compiess  dm  received  sign.,  into  a  compadble  dynamic  range.  Tmdidonally,  airborne 
surveillance  RADARs  have  employed  an  analog  feedback  loop  configuration  ,o  achieve  dm  required 
au,omadc  gain  connol  (AGC)  tmtion.  A  .ypical  implenmnmdon  of  dm  analog  AGO  using  dm  well  known 
ctosed  loop  negative  feedback  configurodon  is  shown  In  fig  I,  when,  Kol)  repreSe„B  ^  fi|Kr 

chatactensuc  of  dm  feedback  am.  The  voltage  condolled  variable  gain  deviee  in  fig.  1  can  be  edher  an 

amplifier  or  adenniato,  or  combinadon  of  bod,  as  determined  by  dm  requimmems  of  die  syslem  umler 
consideration. 


Figure  1:  Typical  analog  AGC  loop  implementation. 

In  this  repon  we  consider  ,he  use  of  a  more  elabonne  and  veisalile  scheme  of  AGC  mad.  possible  by 
ud  rang  the  digital  processing  capabilides  of  the  newly  developed  digital  scan  convener.  To  distinguish 
bmween  die  two  approaches,  the  uadidonal  analog  aulomadc  gain  condol  will  be  subsequently  tefeired  u, 
as  Analog  Automatic  Gain  Condo,  (AAGC)  am,  die  digital  approach  as  Digila,  Automadc  Gain  Condol 
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(DAGC).  The  DAGC  approach  forms  an  estimate  of  the  required  gain  for  each  point  by  building  an 
average  clutter  map  of  the  region  under  consideration  over  the  course  of  multiple  scans.  This  clutter  map 
can  be  digitally  manipulated  to  calculate  the  required  gain  for  returns  from  each  geographic  location.  The 
appropriate  control  voltage  is  then  applied  to  the  controlling  amplifier/attenuator  to  achieve  the  required 
gain.  The  proposed  configuration  of  the  DAGC  is  illustrated  in  figure  2. 

The  report  is  divided  into  4  sections.  Section  2  summarizes  the  theory  describing  the  operation  of  a 
traditional  AGC  system  and  develops  the  control  theory  for  the  proposed  DAGC  system.  In  section  3 
simulations  are  run  to  clarify  the  performance  of  the  AAGC  and  DAGC  systems.  Section  4  summarizes  the 
report  and  makes  recommendations  for  further  work. 

1.1  Receiver  Specifications 

For  the  purposes  of  this  report  simplified  AGC  configurations  are  assumed  in  which  only  one  controlling 
element  is  used  (i.e.  a  voltage  controlled  attenuator).  In  actual  systems,  it  often  proves  necessary  to  place 
gain  control  elements  throughout  the  input  chain  and  frequently  successive  control  loops  are  employed. 
Nevertheless,  the  following  theory  and  results  are  still  applicable  and  it  is  relatively  straightforward  to 
adapt  the  analysis  to  these  more  complex  configurations. 

Table  1.1  summarizes  the  nominal  receiver  specification  assumed  in  this  report.  Where  possible 
specifications  have  been  derived  from  the  AN/APS-506  surveillance  RADAR  to  ensure  a  realistic  model 
has  been  used.  Fixed  amplification  and  attenuation  stages  are  assumed  inserted  throughout  the  system  to 
meet  minimum  and  maximum  signal  levels  as  applicable.  The  controlling  element  is  assumed  to  be  three 
cascaded  Watkins  Johnson  G2  voltage  controlled  attenuators.  The  choice  of  an  attenuator  or  amplifier  will 
obviously  be  system  dependent  but  the  theory  and  results  presented  in  this  report  are  applicable  to  either 
configuration. 

Table  1.1:  RADAR  Simulation  Parameters 


Beamwidth 

2.2° 

Pulse  Repetition  Frequency  (PRF) 

1987  (pulses/sec) 

Scanning  Rate 

300  rpm 

Range  Resolution 

<  lm 

Noise  Temperature 

45  OK 

A  to  D  quantization  level 

0.0039V 

A  to  D  bits  &  range 

8  bits  / 1  volt 

Antenna  Aperture 

1.0m  *  0.5m 

Squint  angle 

15° 

Range 

5000  m 

Swath  width 

4800  m 

Peak  Transmitted  Power 

500  kW 

Noise  Bandwidth 

500  MHz 

Front  end  gain  (pre-AGC) 

100  dB 
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2.0  Automatic  Gain  Control  Theory 


2.1  Analog  Automatic  Gain  Control  (AAGC) 


The  theory  describing  the  performance  and  operation  of  the  traditional  AGC  loop  has  been  well  developed 
[1]  and  only  a  brief  summary  is  given  below.  Figure  3  illustrates  the  AAGC  configuration  assumed  in  this 
report.  The  defining  equations  for  the  AAGC  loop  are  easily  determined  by  inspection  of  figure  3  and  are 
given  by 

c  =  M(V)*r  2.1 

and 

V  =  C  *  fix  *  P{(0)  ,  2.2 


where  r  and  c  represent  the  amplitude  of  the  input  and  output  signals  respectively,  P]  is  the  gain  in  the 
feedback  arm  with  the  envelope  detector  assumed  to  be  operating  in  a  linear  detection  mode,  and  P((o) 
represents  the  filter  characteristic  of  the  feedback  arm  normalized  to  one  at  DC.  M(V)  represents  the 
voltage  amplification  factor  and  is  controlled  via  the  feedback  arm  voltage.  Differentiating  2.1  with  respect 
to  r,  re-arranging  and  substituting  in  2.1  and  the  derivative  of  equation  2.2  with  respect  to  M,  we  obtain 


dc 

c 


1  dr 


1-//1 0(6)) 


c  dM 
M  dV 


r 


2.3 


Linearizing  about  the  DC  operating  voltage,  Vdo  we  obtain  the  frequency  dependent  modulation  transfer 
characteristic 


dc 

c 


1  dr 


1 


Vdc  dM.Dc 
M  dV 


2.4 
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where -  js  the  derivative  evaluated  at  Vdc- 

dV 


The  subtractor  circuit  on  the  attenuator  control  line  in  figure  3  converts  the  attenuator’s  control 
characteristic  to  a  positive  function  i.e.  attenuation  increases  with  increasing  voltage.  The  combined 
control  characteristic  produced  by  the  G2  attenuator  chain  and  subtractor  is  plotted  in  figure  4. 


Three  G2  Voltage 
Controlled  Attenuators 
in  Series 


Figure  3:  AGC  configuration  assumed  in  report 
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Table  2.1:  Regulation  Characteristics 


Input  Power  (W) 

Output  Power  (W) 

Dynamic  range  between 
noise  and  land  clutter  on 
output 

JU  1=10 

Noise 

2.5*10'2 

7.8*10'5 

Sea  Clutter 

6.5*10‘' 

1.8*10 3 

Land  Clutter 

6.5*10' 

2.4*1  O'2 

25  dB 

H  ,=150 

Noise 

2.5*10’2 

7.8*  10'5 
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The  static  regulation  curve  of  steady  state  input  versus  steady  state  output  can  now  be  plotted  by  calculating 
the  voltage  for  a  given  output  signal  c  using  the  known  value  of  pi,  and  then  reading  the  corresponding  gain 
from  figure  4  to  calculate  r.  (3(0)  is  assumed  equal  to  1.  Figure  5  displays  the  regulation  curves  for  a 
variety  of  pi  values.  Table  2.1  summarizes  the  input  and  output  power  associated  with  receiver  noise,  land 
clutter  and  sea  clutter  for  the  three  different  feedback  arm  gains  shown  in  figure  5.  All  noise  is  assumed  to 
be  Rayleigh  distributed  with  the  mean  scattering  cross-section  density  assumed  to  be  0.01  and  0.0001 
m2/m2  for  land  and  sea  clutter,  respectively.  The  input  powers  are  calculated  using  the  parameters  listed  in 
table  1.1  and  are  also  used  for  the  simulations  carried  out  in  section  3.0.  In  all  cases  the  impedance  has 
been  assumed  to  be  1  ohm.  The  choice  of  front  end  gain  is  arbitrary  (i.e.  gain  prior  to  AGC  loop),  the 
important  factor  is  where  on  the  static  regulation  curve  the  AAGC  is  operating.  As  can  be  seen  in  figure  5, 
a  decrease  in  the  front  end  gain  can  be  compensated  for  by  increasing  the  feedback  arm  gain,  pb  which 
results  in  the  static  regulation  curve  being  shifted  along  the  slope  defined  by  the  unregulated  portion  of  the 
curve.  By  judicious  choice  of  pi,  the  same  relative  regulation  between  input  and  output  can  be  achieved 
although  the  absolute  value  of  the  input  and  output  will  be  altered  in  proportion  to  the  change  in  the  front 
end  gain  and  thus  can  be  adjusted  to  meet  the  power  level  requirements  of  the  system  components.  All 
outputs  are  normalized  so  that  the  average  voltage  output  due  to  receiver  noise  alone  is  equal  to  two  A/D 
quantization  units  or  0.0078  V.  Using  table  2.1  we  choose  a  pi  value  of  150  as  providing  a  reasonable  level 
of  regulation  i.e.  the  34  dB  dynamic  range  of  input  between  receiver  noise  and  land  clutter  is  compressed  to 
5  dB  on  output.  This  value  will  be  used  for  subsequent  calculations  involving  the  AAGC  performance. 
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Output  Power  (dBm)  3  Attenuation  (dB) 


The  frequency  behaviour  of  the  AAGC  can  now  be  calculated  using  equation  2.4.  In  this  example,  the  low 
pass  filtering  action  in  the  feedback  arm  is  assumed  to  be  accomplished  using  a  range  gated  filter  (RGF) 
bank  of  1.0  pF  capacitors  attached  in  parallel  to  the  feedback  arm  as  shown  in  figure  3.  The  equivalent 
circuit  is  that  of  a  simple  RC  lowpass  filter.  This  configuration  is  identical  to  the  RGF  currently  employed 
in  the  AN/APS-506  RADAR.  Consistent  with  the  AN/APS-506,  each  capacitor  is  assumed  gated  on  for  a 
range  interval  of  3.8  ps,  and  switched  to  a  hold  configuration  (i.e.  input  line  floating)  during  the  ‘off 
periods. 


Figure  6:  Sampling  action  of  range  gated  filter  banks.  Dotted  line  shows  impulse  response  of  RC 
filter  in  time  and  dashed  lines  show  how  impulse  response  is  distributed  in  time  by  range  gating 
action.  Solid  line  is  a  sample  sinusoid  input  signal.  Current  output  of  RGF  (TIME  =  0)  would  equal 
integrated  sum  of  dashed  filter  response  multiplied  by  sample  signal  from  TIME=0  through  to  (- 
Infinity). 

Figure  6  demonstrates  the  applicable  sampling  action  where  the  time  domain  impulse  response  of  the  RC 
filter  has  been  distributed  across  the  discrete  time  samples  by  the  holding  action  of  the  RGF.  Figure  7 
displays  the  lowpass  behaviour  of  this  range  gated  filter  calculated  assuming  that  each  ‘on’  period  can  be 
modeled  as  an  impulse  function.  This  assumption  is  valid  as  long  as  the  wavelength  is  much  longer  than 
the  length  of  the  range-sampling  period.  Errors  of  less  than  2.5%  are  introduced  for  wavelengths  8  times  as 
long  or  greater  than  the  sampling  window.  The  corresponding  AAGC  transfer  functions  are  also  shown  for 
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a  variety  of  input  signal  amplitudes  assuming  Pi=150,  Because  an  undelayed  AAGC  is  assumed  (i.e.  no 
reference  voltage)  lower  input  power  signals  are  more  poorly  regulated  than  higher  power  signals.  The 
overall  frequency  behavior  of  the  AAGC  system  is  that  of  a  high  pass  filter,  as  expected. 


Figure  7:  Bode  plot  of  range  gated  filter  bank  and  AAGC  loop  with  feedback  arm  gain,  pi,  equal  to 
150.  Curves  are  displayed  for  average  AAGC  input  signal  power  levels  after  input  amplification  and 
are  as  follows:  receiver  noise  =  .025  W,  sea  clutter  =  0.65  W  and  land  clutter  =  65  W. 

Clearly  the  performance  of  the  AAGC  is  limited  by  the  fixed  nature  of  the  RGF,  i.e.  the  capacitors  are 
fixed,  and  also  the  relatively  slow  rolloff  of  the  AAGC  frequency  response.  In  particular,  the  long  memory 
of  the  AAGC  will  degrade  the  performance  of  the  gain  control  in  viewing  configurations  where  the 
azimuthal  movement  of  the  beam  is  from  a  strong  clutter  environment  such  as  land  clutter  to  a  weaker  sea 
clutter  environment. 

2.2  Digital  Automatic  Gain  Control  (DAGC) 

Figure  2  displays  a  block  diagram  describing  the  proposed  operation  of  the  DAGC.  RADAR  returns 
entering  the  DAGC  are  pulse  compressed  and  mapped  to  their  polar  coordinates  in  memory.  This  polar 
coordinate  data  is  then  scan  converted  to  a  raster  format  to  form  a  clutter  map  of  the  returns.  A  typical  scan 
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conversion  algorithm  is  described  in  detail  in  DREO  Report  No.  1113  [2]  and  involves  the  filtering  of  the 
data  to  prevent  aliasing  followed  by  interpolation  to  the  raster  grid.  The  raster  data  is  now  available  for 
manipulation  by  the  processors  in  the  DAGC  feedback  arm.  The  first  stage  of  manipulation  allows  the 
addition  of  multiple  scans  to  improve  the  signal  to  noise  ratio  of  the  retrieved  clutter  map.  The  averaged 
raster  data  is  then  fast  fourier  transformed  (FFT)  from  the  spatial  domain  to  the  frequency  domain  where  it 
can  be  digitally  filtered  as  desired  by  multiplication  with  a  2-D  filter  function.  The  filtered  clutter  map 
values  are  now  compared  with  the  desired  threshold  voltage  to  determine  the  gain  correction  required.  In 
this  report  the  threshold  is  set  equal  to  two  quantization  levels  of  the  A  to  D  converter.  This  represents  a 
compromise  between  maximizing  the  dynamic  range  available  for  detection  while  ensuring  quantization 
loss  is  kept  to  a  reasonable  level.  As  a  final  step  the  generated  gain  control  vector  is  decompressed  for 
application  to  the  adjustable  attenuator  which  is  located  prior  to  the  pulse  compression  unit  in  the  input 
chain. 


2.2.1  Transfer  Function  and  Stability  of  DAGC  without  Spatial  Filtering 

Since  the  DAGC  is  a  digital  system  sampling  the  input  signal  once  per  time  step,  the  transfer  function  of 
the  loop  is  most  easily  evaluated  using  the  z  transform.  This  can  be  determined  by  taking  the  Laplace 
transform  of  the  periodic  discrete  signal  and  replacing  the  factor  esT  with  z,  where  T  is  the  sampling  interval 
of  the  signal. 


A  more  detailed  description  of  z-transforms  and  their  applications  can  be  found  in  [3].  It  can  be  show  that 
the  transfer  function  of  a  delay  element  of  one  time  step  is  z'1.  In  order  to  make  the  analysis  of  the  system 
more  amenable  to  standard  feedback  theory,  it  is  convenient  to  work  with  the  log  of  the  signals  and  convert 
all  operations  to  their  log  equivalents.  Using  the  above  definition,  the  flow  chart  of  a  simplified  system 
without  spatial  filtering  is  shown  in  figure  8.  This  simplified  view  is  equivalent  to  modeling  the  DAGC 
properties  for  one  geographic  location  independent  of  the  points  surrounding  it. 


Figure  8:  Flow  chart  of  simplified  DAGC  loop  without  spatial  filtering.  N  equals  the  number  of 
temporal  samples  averaged  together  and  REF  is  the  desired  reference  level. 
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Each  state  variable  is  identified  by  a  subscript  number  after  the  x.  Intermediate  variables  are  also  shown  in 
the  figure,  these  variables  are  identified  by  a  letter  subscript.  These  variables  are  redundant  and  are  not 
required  for  the  matrix  solution,  they  are  introduced  simply  to  facilitate  the  discussion. 

In  a  real  system  the  averaging  process  is  applied  to  a  non-log  signal,  i.e. 

(s  +  t/l)  +  Cs  +  t/2)  +  (.s  +  d3)  _  dl  +  d2  +  d3 

y  =  — — - — - =  S  H - . 

3  3 

where  we  have  temporarily  set  REF  =  0  for  simplicity.  Assuming  that  the  fluctuations,  d;  are  equal  and 
represent  the  independent  standard  deviations  of  probability  distributions  governing  the  noise  they  are  then 
added  in  quadrature  giving 


Taking  the  log  we  obtain 

log(*J  =  logO  +  -J=0.  2/7 

In  the  log  converted  flow  chart  we  have  approximated  this  averaging  process  by  averaging  the  log  signal 
directly  to  give  LOGxa  which  can  be  shown  to  be  roughly  equivalent  to  logfxj  provided  the  fluctuations  in 
d  are  small  compared  to  the  average  signal  level  s.  Averaging  in  the  log  system  gives  us 

_  [logO  +  ^1)  +  logO  +  +  logO  + 1/3)]  2  S 

L/L/IJa  “  * 


But 


d(l0g  5  therefore  if  d;  is  small  compared  to  s  we  can  linearize  about  s  to  obtain 

9  s  s 


LOGxa  ~ 


log(s)  +  —  +  log(s)  +  ^  +  logO?)  +  ^  dl  d2  d3 

- £ - —1 - =  log(s)+_+_  +  _. 


2.9 


As  above  the  fluctuations  d;  add  in  quadrature  giving 
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2.10 


d  d 

LOGxa  =  log(s)  +  — j=  =  log(s  +  -j=)  =  log(.xa ) . 

The  feedback  loop  in  figure  8  can  now  be  analyzed  using  standard  techniques.  The  families  of  feedback 
loops  present  in  the  flow  chart  are  identified  in  figure  8  as  LI,  L2  and  L3  and  are  listed  below. 


z(N-(N-0))  *_L*(z-1)*(_A) 

N 

1 

1 

_ i 

z(W))*J_*(z-l)  +  (_A) 
N 

— 

-z~2A 

N 

-Z'NA 

• 

N 

1 

Nl 

l 

a* 

_  _ 

N  . 

L2  =  [z_1  J  2.12 


u  = 


(-z-'z~(N-2)  L)*(-z-V("-3)  1) 

N  N  N 

--(w-3)  _L\*  */'_7-u-(amV)  _L\ 

1  n}'"k  nj 


N 


2.13 


where  A  is  a  variable  introduced  to  represent  the  accuracy  to  which  the  gain  can  be  adjusted  in  an  actual 
system,  i.e.  A=1  would  represent  the  ability  to  apply  perfect  corrections.  Using  Mason’s  rule  [4]  the 
transfer  function  from  r  to  c  is 


T(r-  >  c)  = 


n  aq 

A 


2.14 


where 

A  =  1  -  2  A,  +  L2  +  X  L3 ,  Ak  =  A  +  X  L, ,  and  P  =  1  for  forward  path. 

If  the  transfer  function  is  to  be  stable  with  time  then  the  roots  of  the  denominator  must  all  be  less  than  or 
equal  to  zero.  Table  2.1  shows  the  stable  range  of  A  values  for  different  values  of  N.  Since  the  voltage 
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control  curve  for  the  controlling  attenuator  is  well  known,  the  variable  A  will  typically  be  quite  close  to  1 
and  the  D  AGC  will  be  stable  for  all  values  of  N. 


Table  2.1  Stable  range  of  A  for  varying  N. 

N _ Range  of  A  for  Stable  Operation 

2  0-3 

5  0-6 

10  0-11 

2.2.2  Noise  Performance  of  DAGC 

The  output,  Cj,  from  the  DAGC  after  j  steps  can  be  written  as 

Cj  =  t j  -  Ar V_,  +  A(A - 1  )r._2  -  A(A  - 1)2  r;_3  +  A(A  - 1)3 r,_4 

+  ....(-l)7-1  A(A  -  l)j~2  rx  2.15 

=0+g(-irA(A-ir‘0, 


where  q  represents  the  input  at  each  step  j. 

To  calculate  the  contribution  from  the  independent  input  noise  at  each  step  to  the  output  signal  we  separate 
ij  into  its  signal  and  noise  components,  rj=s+dj.  As  above,  s  should  be  a  constant  value  at  any  one  location 
and  dj  equals  the  standard  deviation  of  the  input  noise  about  this  value.  Linearizing  around  s  as  shown 
above  and  adding  the  noise  components  quadratically  we  obtain 


d'out 


d 


■  +  A‘ 


j  2  ,2 

-^-  +  A2(A-l)2-^- 


+  A  [(A-l)  ] 


2  1 2  ^  j~ 3 


•  +  A2  [(A  - 1)3]2  —^-7—  + ... 


XA2[(A-l)m-1]2 

m= 1 


2.16 


This  series  quickly  converges  for  values  of  A  within  the  stability  limits  of  the  DAGC  system.  It  should 
again  be  emphasized  that  dout  is  the  log  value  of  the  deviation  due  to  noise,  i.e.  if  the  log  of  the  constant 
portion  of  the  output  signal  is  sout  then  the  log  output  will  have  a  one  sigma  variation  about  sout  such  that 
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c=sout  +/-  dout.  Typically  the  noise  figure  is  constant  so  that  dj=d.  If  the  attenuator  voltage  characteristic  is 
also  perfectly  known  then  A=1  and  equation  2.16  reduces  to 


2.17 


Figure  9  shows  a  plot  of  the  output  noise  enhancement  factor  (noise  in/  noise  out)  from  the  DAGC  versus 
A.  Figure  9  was  produced  by  iteratively  adding  numerically  generated  guassian  noise.  Since  the  log 
approximation  is  not  used  the  results  in  figure  9  will  differ  slightly  from  that  those  generated  using  equation 
2.16.  Also  shown  are  two  curves  for  the  case  in  which  3  and  10  scans  are  averaged.  The  improvement  in 
noise  performance  that  can  be  achieved  by  scan  averaging  is  demonstrated  for  all  values  of  A.  The 
introduction  of  scan  averaging  also  increases  the  response  time  of  the  feedback  loop.  Figure  10 
demonstrates  the  response  of  a  DAGC  system  to  a  step  input  with  A  ranging  from  0.75  through  to  1.25  for 
the  case  of  1,  3  and  10  averaged  scans. 


0.75  0.85  0.95  1.05  1.15  1.25 

A 

Figure  9:  Noise  enhancement  factor  versus  accuracy  with  which  correction  can  be  applied.  A=1 
represents  perfect  application  of  desired  gain  correction  while  A<1  represents  an  underestimate  and 
A>1  an  overestimate  of  the  required  correction.  Curves  are  shown  for  cases  in  which  1, 3  or  10  past 
scans  are  averaged  to  form  the  gain  correction  estimate. 
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2.2.3  Frequency  Responses  of  DAGC 


The  frequency  response  of  the  transfer  function  can  be  determined  using  the  definition  of  z  and  the 
transformation 


i  +  Mv 

i  -  Mv 


2.18 


where 


,Q)T 

cow  =  tan(— — ) 


and  to  is  the  angular  frequency  of  the  signal. 


Figure  11:  Bode  plot  of  DAGC  response  for  N=l,  3  and  10  scan  averaging,  x-axis  is  actually  w*T,  to 
get  frequency  divide  the  x-axis  values  by  the  sampling  period,  T. 
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The  frequency  responses  of  the  loop  for  different  values  of  N  are  shown  in  fig  11.  The  curves  in  fig.  1 1 
show  that  as  the  number  of  scans  averaged  is  increased  the  cutoff  of  the  frequency  response  is  pushed  to 
lower  frequencies.  It  should  be  emphasized  that  this  frequency  response  refers  to  the  temporal  behaviour  of 
the  DAGC  for  the  returns  from  one  geographic  location  i.e.  from  scan  to  scan  how  does  the  return  at  one 
point  on  the  clutter  map  vary  with  time.  As  might  be  expected  the  frequency  response  displays  a 
periodicity  with  to  which  could  result  in  aliasing  for  an  insufficiently  sampled  figure. 

2.2.4  Transfer  Function  and  Stability  of  DAGC  with  Spatial  Filtering 

The  introduction  of  spatial  filtering  across  the  imaged  scene  rapidly  increases  the  complexity  and 
computational  load  required  to  solve  for  the  transfer  function.  A  state  variable  technique  is  used  below  to 
solve  the  problem.  Figure  12  illustrates  the  flow  chart  for  one  arm  of  the  overall  flow  chart  corresponding 
to  the  input  to  output  for  one  location,  i,  on  the  clutter  map.  Similar  flow  charts  will  exist  for  each  location 
on  the  clutter  map.  The  arms  are  coupled  together  through  weighted  averaging  in  the  feedback  section  of 
each  arm.  There  are  N+l.  state  variables  associated  with  each  arm.  Each  variable  is  identified  by  a  number 
after  the  x  and  a  subscript  to  identify  the  arm  or  spatial  location  to  which  the  state  variable  belongs. 
Intermediate  variables  are  also  shown  in  the  figure,  these  variables  are  identified  by  a  letter  subscript. 
These  variables  are  redundant  and  are  not  required  for  the  matrix  solution,  they  are  introduced  simply  to 
facilitate  the  discussion.  The  coupling  to  other  arms  occurs  at  point  xbj  where  the  scan  averaged  values  are 
fed  to  the  other  arms  of  the  overall  flow  chart.  Likewise,  variable  x^  corresponds  to  the  value  formed  from 
the  weighted  addition  of  xb>i  with  the  xb  values  from  other  arms.  The  choice  of  weighting  function  applied 
determines  the  spatial  filtering  applied  to  the  feedback  signal.  In  the  actual  systems  this  step  would 
probably  be  accomplished  by  the  computationally  more  efficient  method  of  performing  a  2D  FFT  of  the 
map  of  xb  variables  and  then  applying  a  2D  filter  in  the  frequency  domain. 

The  general  form  of  the  dynamic  equations  for  the  system  is 

x[(j+l)T]  =  Ax(jT)  +  Bu(jT)  2-19 

c(jT)  =  Dx(jT)  +  Eu(jT)  2,20 

where  x  is  the  matrix  of  state  variables,  u  is  the  matrix  of  inputs,  c  is  the  matrix  of  outputs  and  A,  B,  D  and 
E  are  defined  below.  In  the  above  j  represents  the  time  step  index.  In  the  subsequent  equations  the 
subscript  i  will  be  identified  with  the  branch  of  the  output  variable  under  consideration  while  k  will  be 
associated  with  the  branch  of  the  variable  acting  as  an  input  to  the  system.  N  represents  the  number  of 
scans  averaged  on  any  one  branch  and  m  represents  the  number  of  branches  or  inputs  into  (and  in  this 
analysis,  outputs  from)  the  system.  REF  represents  the  desired  reference  level. 
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Figure  12:  Flow  chart  of  simplified  DAGC  loop  with  spatial  filtering.  See  text  for  description  of 
variables  and  subscripts.  See  figure  8  for  loop  definitions. 

Each  of  the  matrices  in  equations  2.19  and  2.20  can  be  determined  by  inspection  of  the  flowchart  in  fig.  12 
and  are  as  follows,  x  is  given  by 


where  xk  is 


xlk  through  xN+lk  represent  the  state  variables  for  branch  k.  A  has  the  dimensions  m*(N+l)  by  m*(N+l) 
and  can  be  formed  by  examination  of  figure  12  giving 


A  = 


al,l 

ah2 

*  a\ ,m 

a2,l 

U2,2 

a2,m 

am,l 

am,2  • 

*  am,m 

where  ai)k  has  the  dimensions  N+l  by  N+l  and  is  given  by 
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The  subscript  on  the  w  identifies  the  state  variable  it  is  associated  with.  aoff  also  has  the  dimensions  N+l  by 
N+l  and  is  given  by 
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and 

aoff  =  0  if  i^k. 

In  the  above  it  has  been  assumed  that  all  branches  use  the  same  scan  to  scan  averaging  coefficients  wv  In 
addition  both  (0  and  w  coefficients  are  assumed  to  have  been  normalized  by  their  respective  factors  of  l/o>r 
and  l/wT  as  shown  in  figure  12. 

B  has  the  dimensions  m*(N+l)  by  m+1  and  is  given  by 
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where  k  refers  to  the  column  of  matrix  bi,  8itk=l  if  i=k  and  6i>k=0  if  i^k.  u  is  an  m+1  by  1  matrix  and  is 
given  by 


h 

r2 
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D  is  given  by 


D  =  [D,  Z>2  ..  Dm] 


where  Dkhas  the  dimensions  m  by  N+ 1  and  is  given  by 
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where  the  subscript  i  refers  to  the  row  of  the  matrix,  8i  t  =1  if  i=k  and  5i  k  =0  if  i*k.  The  matrix  E  has  the 
dimensions  m  by  m+1  and  is  given  by 


1  0  ..  0  0 

0  1  ..  0  0 

0  0  ..  1  0 


It  can  be  shown  [3]  that  the  transfer  function  of  the  system,  G=c/r  is  given  by 
G(z)=[D(zI-A  )-‘B+E]. 


2.21 


Once  the  transfer  function  of  the  system  has  been  determined  the  stability  of  the  system  can  be  affirmed  in 
the  usual  manner  by  checking  if  the  absolute  roots  of  the  denominator  in  equation  2.21  are  less  than  one  i.e. 

|z|<  1. 
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3.0  Simulations 


In  order  to  illustrate  the  differences  between  an  AAGC  and  DAGC  system,  a  series  of  simulations  were 
performed.  As  stated  earlier  it  is  anticipated  that  the  main  advantage  of  the  DAGC  system  will  occur  along 
regions  such  as  shorelines  where  the  AAGC  can  suffer  residual  affects  from  the  large  scattering  cross 
section  of  the  land  clutter  as  the  beam  sweeps  from  the  land  onto  the  ocean  surface.  Figures  13  and  14 
illustrate  the  two  viewing  configurations  simulated  in  this  section.  The  first  configuration  is  chosen  to 
simulate  a  sector  viewing  operation  where  the  beam  sweeps  from  sea  onto  land  and  back  onto  sea  again. 
The  second  configuration  simulates  the  effect  of  the  beam  sweeping  parallel  to  the  shore  with  the  range 
bins  corresponding  to  the  target  range  remaining  equidistant  from  the  land-sea  transition.  The  average 
range  in  both  cases  is  assumed  to  be  5  km  and  the  change  of  azimuthal  sample  spacing  with  range  is 
insignificant  (i.e.  has  little  effect  on  the  results)  over  the  swath  width  under  consideration  and  is  assumed 
constant  for  these  simulations.  To  further  simply  the  simulations  other  complicating  factors  such  as  pulse 
compression,  peak  detection,  beam  width  etc.  are  ignored.  Range  dependence  of  the  returns  across  the 
swath  width  is  also  ignored. 

All  clutter  in  these  simulations  is  assumed  to  be  Rayleigh  distributed  with  the  mean  scattering  cross  section 
varied  as  necessary  to  represent  land  versus  sea  clutter.  Land  clutter  is  assumed  to  have  a  scattering  cross 
section  density  of  0.01  m2/m2  while  sea  clutter  is  assumed  equal  to  0.0001m2/m2.  The  discrete  targets  in 
the  scene  are  assumed  to  have  a  cross  section  of  1  m2.  To  reduce  the  computational  load,  range  samples  are 
generated  every  75  m  i.e.  64  range  intervals  over  the  4800  m  wide  swath.  However  a  range  interval  length 
of  less  than  one  meter  is  still  assumed  during  the  calculation  of  the  clutter  viewing  area  and  the 
corresponding  clutter  signal  strengths  to  preserve  the  clutter  to  signal  ratios  which  would  be  expected  in  an 
actual  high  resolution  RADAR.  For  each  case  the  equivalent  of  200  scans  were  generated,  with  the  random 
noise  generator  re-seeded  for  each  scan.  These  scans  are  consecutively  processed  as  per  the  appropriate 
algorithm  to  simulate  the  AAGC  and  DAGC  systems.  The  applicable  algorithms  for  the  AAGC  and  DAGC 
simulators  are  shown  in  the  block  diagrams  of  figures  15  and  16  respectively.  The  final  output  is  generated 
using  the  simple  recursive  formula 

Output  Image  =  a*new  scan  +  (l-oc)*old  scan  3.1 

with  a-0.025. 

For  all  cases  an  appropriate  amplification  has  been  applied  to  the  AAGC  output  to  ensure  that  that  the 
output  due  to  receiver  noise  alone  would  equal  two  quantization  units  of  an  8  bit,  1  volt  range  AtoD  (i.e. 
1/255  V).  The  same  desired  quantization  level  is  automatically  implemented  within  the  threshold 
comparison  step  of  DAGC  processing  where  the  threshold  has  been  specified  at  the  same  value. 
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42  bins  172  42  bins 


Figure  13:  Test  configuration  1  to  simulate  the  effect  of  a  radar  beam  sweeping  from  sea  clutter  onto 
land  and  back  onto  sea  clutter.  Azimuthal  bin  width  is  79m  and  range  bin  width  is  75m. 
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256  bins 


Figure  14:  Test  configuration  2  to  simulate  the  effect  of  radar  fc 
AAGC  sampling  window  of  3.8  ps  straddling  land-sea.  Azimutl 
width  is  75m. 
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3.1  Configuration  1  Simulations  -  RADAR  Sweeping  from  Sea  to  Land  to  Sea 


Figure  17  illustrates  the  output  from  the  AAGC  loop  for  configuration  1  input  data.  To  clarify  the 
discussion  the  return  signal  from  the  land  and  sea  clutter  has  been  divided  into  a  deterministic  portion 
corresponding  to  the  mean  return  expected  from  the  land  and  sea  clutter  and  a  non-deterministic  portion 
corresponding  to  the  fluctuations  about  this  mean.  Figure  17  shows  only  the  deterministic  portions  of  the 
clutter  signals  and  the  target  signals.  The  AAGC  loop  configuration  assumed  was  described  previously  and 
is  illustrated  in  figure  3,  the  capacitance  and  resistance  of  the  RC  filter  components  are  1  pF  and  392  £2 
respectively  with  three  cascaded  G2  attenuators  providing  the  gain  control.  Unless  otherwise  specified  the 
feedback  arm  gain,  p,1(  is  150. 

The  movement  of  the  RADAR  beam  is  from  left  to  right  in  figure  17,  the  lag  in  the  response  of  the  AAGC 
on  the  trailing  edge  of  the  land  to  sea  clutter  transition  is  readily  apparent  with  significant  masking  of  the 
targets  located  at  approximately  17.1  and  17.9  km  (azimuthal  bin  locations  217  and  227).  For  this  example 
the  AAGC  output  does  not  recover  to  the  90%  steady  state  value  until  the  beam  has  moved  a  distance  of 
approximately  2.4  km  or  30  azimuthal  bins  from  the  land  to  sea  transition. 

Examination  of  the  leading  edge  for  the  sea  to  land  transition  reveals  a  similar  lagging  response  as 
evidenced  by  the  large  spike  starting  at  3.3  km  (azimuthal  bin  43).  However  since  this  lag  occurs  over  the 
land  where  we  are  not  searching  for  targets  it  does  not  degrade  the  detection  performance,  in  fact  due  to  the 
non-anticipatory  nature  of  the  AAGC  and  its  simple  RC  filter  we  achieve  single  bin  azimuthal  resolution  of 
targets  on  the  leading  edge. 

The  difficulty  of  identifying  targets  in  the  trailing  zone  of  the  land  to  sea  transition  is  illustrated  in  figure 
18  where  the  non-deterministic  portion  of  the  sea  and  land  clutter  responses  has  been  added  to  the  output 
signal.  It  is  clear  in  figure  18  that  the  non-deterministic  portion  masks  the  already  weak  target  responses 
seen  in  figure  17.  The  plotted  output  is  derived  using  an  alpha  of  0.025  in  equation  3. 1. 

Figure  19  displays  the  output  from  the  DAGC  loop.  It  should  be  noted  that  the  DAGC  process  is  actually  a 
two  dimensional  one  as  the  gain  matrix  correction  applied  is  generated  by  performing  a  two  dimensional 
filtering  operation  on  the  clutter  map  built  up  over  the  previous  scans.  In  the  following  examples  we  have 
assumed  a  simple  2-D  circular  box  filter  with  a  radius  of  200  m.  For  clarity  we  have  extracted  one 
azimuthal  vector  at  the  target  range  from  the  overall  output  matrix.  The  most  obvious  characteristic  of  the 
DAGC  response  is  the  symmetrical  response  of  output  between  the  two  land-sea  interfaces.  This  is 
expected  due  to  the  spatial  filtering  operation  on  the  clutter  map.  The  use  of  a  200m  box  filter  results  in  a 
clear  identification  of  targets  on  the  leading  and  trailing  edges  of  the  land-sea  transitions,  including  the 
targets  located  225m  from  shore.  This  is  in 
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continue  loop 
for  j=l:n 
where  n  is 
total  number 
of  time  steps 


continue  loop 
for  k=l:N 
where  N  is 
total  number 
of  scans 


Figure  15:  Flow  chart  for  algorithm  corresponding  to  AAGC  processing  module. 
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continue  loop 
for  j=l:N 
where  N  is 
total  number 
of  scans 


Figure  16:  Flow  chart  for  algorithm  corresponding  to  DAGC  processing  module. 
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Figure  17:  Output  from  AAGC  for  configuration  1  deterministic  input  signal.  Arrows  identify 
targets. 


Figure  18:  Output  from  AAGC  for  configuration  1  with  non-deterministic  input  signal  added. 
Arrows  identify  targets. 
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0.07 


Figure  19:  Plot  along  azimuthal  vector  containing  targets  of  DAGC  output  for  configuration  1  with 
deterministic  input  data.  Arrows  indicate  targets. 

contrast  to  thenon-symmetric  response  of  the  AAGC  where  clear  identification  of  the  targets  right  up  to  the 
land-sea  transition  on  the  leading  edge  is  achieved  at  the  cost  of  significant  masking  on  the  trailing  edge.  It 
is  also  apparent  that  targets  closer  to  the  land-sea  transitions,  i.e.  within  the  filter's  200  m  radius,  would 
suffer  from  masking  under  the  DAGC  systems.  Roughly  speaking,  spatial  filtering  on  the  clutter  map 
under  the  DAGC  scheme  rather  than  then  the  real  time  approach  of  the  AAGC  allows  the  response  lag  to  be 
symmetrically  distributed  across  the  land-sea  boundary  thereby  cutting  in  half  the  extent  of  the  shadow 
region  on  the  trailing  edge  if  an  identical  filter  width  is  assumed  in  both  approaches.  Of  course,  as 
mentioned  earlier,  this  comes  at  the  expense  of  the  addition  of  a  shadow  zone  on  the  leading  edge  of  the 
land-sea  transition.  In  addition,  it  is  unlikely  that  the  DAGC  and  AAGC  will  have  identical  filter  responses 
although  due  to  the  digital  nature  of  the  DAGC  system  it  is  possible  to  choose  and  implement  any  desired 
filter  response. 
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The  non -symmetrical  AAGC  and  symmetrical  DAGC  behaviour  was  anticipated  and  was  the  primary 
motivation  for  investigating  the  use  of  the  DAGC.  The  resolution  of  the  AAGC  system  can  be  improved 
by  decreasing  the  capacitance  of  the  RC  filter.  Figure  20  illustrates  the  effect  of  reducing  the  capacitance 
by  two  orders  of  magnitude  to  0.01  pF.  Practically,  the  manipulation  of  the  capacitance  of  the  range  gated 
filter  banks  is  very  difficult  under  actual  operation  and  it  is  not  possible  to  adjust  the  capacitance  to  match 
the  varying  pulse  repetition  and  scanning  rates  used  in  the  different  modes  of  operation.  This  typically 
results  in  the  choice  of  a  compromise  value  and  the  resulting  non-optimal  performance  of  all  modes.  In 
contrast,  manipulation  of  the  DAGC  parameters  is  trivial  as  the  filtering  operation  is  performed  digitally 
and  operating  parameters  are  easily  adjusted  within  the  software. 
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Figure  20:  Output  from  AAGC  for  configuration  1  deterministic  input  signal  with  capacitance 
reduced  by  two  orders  of  magnitude  in  RC  filter.  Arrows  indicate  targets. 

Figure  21  shows  the  output  with  the  non-deterministic  portion  of  the  clutter  signal  added  using  an  alpha  of 
0.025  in  equation  3.1  above.  Comparing  figures  18  and  21  the  performance  improvement  on  the  trailing 
edge  and  small  degradation  on  the  leading  edge  under  DAGC  processing  is  readily  apparent. 

Another  desirable  feature  of  the  DAGC  approach  is  the  excellent  regulation  of  the  return  signal  across  the 
full  field  of  view.  Other  than  small  transition  regions  where  the  box  filter  straddles  the  sea  land  interface, 
the  output  is  perfectly  flat  (ignoring  the  non-deterministic  portion  of  the  clutter)  and  exactly  equal  to  the 
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desired  quantization  level,  which  is  that  due  to  the  intrinsic  receiver  noise  temperature  alone.  By  contrast 
the  output  from  the  AAGC  loop  is  imperfectly  regulated,  and  while  the  dynamic  range  between  the  sea 
clutter  and  land  clutter  is  greatly  reduced  there  is  still  a  1.2  db  difference  and  in  both  cases  the  output  is 
above  that  due  to  receiver  noise  alone.  While  the  regulation  response  of  the  AAGC  can  be  improved  by 
increasing  the  gain  in  the  feedback  arm  (or  alternately  by  increasing  the  input  signal  strength  into  the  loop 
see  table  2.1),  this  requires  the  loop  to  be  operated  at  a  higher  attenuation  level.  This  may  exceed  the 
maximum  dynamic  range  of  the  voltage  controlled  attenuators  (or  equivalently  amplifiers)  and  also  require 
adjustment  to  the  other  gains  in  the  input  chain.  This  type  of  adjustment  proves  difficult  in  actual  systems 
and  introduces  the  possibility  of  operator  error. 


Figure  21:  Plot  along  azimuthal  vector  containing  targets  of  DAGC  output  for  configuration  1  input 
data  with  non-deterministic  portion  of  clutter  signal  added.  Arrows  indicate  targets. 


It  should  be  noted  that  the  improved  resolution  of  the  DAGC  approach  implemented  in  this  report  comes  at 
the  expense  of  reduced  detection  probability  in  regions  where  the  fine  resolution  is  not  required.  This 
occurs  due  to  the  small  spatial  extent  of  the  ‘box’  filter  used  to  extract  the  gain  correction  signal  from  the 
clutter  map.  The  smaller  the  spatial  dimensions  of  the  filter  used  the  more  heavily  the  actual  target  will 
contribute  to  the  gain  correction  matrix  in  the  region  around  the  target.  Hence  the  gain  in  this  region  will 
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be  pulled  down  by  the  presence  of  the  target  thereby  reducing  the  final  output  signal  strength  of  the  target 
and  its  detection  probability.  This  is  also  true  for  the  AAGC  system.  If  the  time  constant  of  the  RC  filter  is 
decreased  to  allow  the  response  to  adjust  more  quickly  to  a  land-sea  transition  then  more  of  the  high 
frequency  components  associated  with  the  target  signal  will  be  suppressed  by  the  AAGC  resulting  in 
similar  reductions  in  the  signal  target  strength. 

3.2  Configuration  2  Simulations  -  RADAR  Sweeping  Parallel  to  Land 

Configuration  2  investigates  the  effect  of  the  finite  sampling  window  of  the  range  gated  filtering  employed 
by  the  AAGC  system.  Each  range  gate  filter  is  assumed  to  sample  the  incoming  signal  for  3.8  pis, 
corresponding  to  a  range  length  of  approximately  570  m.  For  the  near  target  located  225  m  from  shore  the 
worst  case  scenario  will  occur  when  the  range  sampling  interval  lies  approximately  half  on  land  and  half  on 
sea  surface  thereby  just  containing  the  target  within  its  limits.  Under  these  conditions  it  is  sufficiently 
accurate  to  assume  that  the  feedback  signal  after  the  RGF  will  equal  the  DC  average  across  the  range 
sampling  interval  as  the  higher  frequency  components  of  the  variation  across  this  window  are  heavily 
filtered  by  the  RGF. 

Figure  22  displays  the  AAGC  outputs  along  a  range  vector  through  targets  located  225  and  975  m  from  the 
shoreline.  The  worst  case  scenario  outlined  above  has  been  assumed.  The  transition  between  range 
sampling  intervals  lying  entirely  over  the  sea  and  the  one  straddling  the  land  sea  interface  is  clearly  visible 
in  the  drop  in  signal  levels  between  the  two  traces.  This  figure  can  be  contrasted  with  figure  23  for  the 
DAGC.  Since  the  range  and  azimuth  bin  intervals  of  the  digital  data  are  similar  in  size,  and  the  filtering 
operation  is  symmetrical  in  the  range  and  azimuth  directions,  the  DAGC  output  of  figure  23  is  essentially 
the  same  as  that  shown  in  figure  21.  Theoretically,  the  azimuthal  resolution  of  the  AAGC  method  could  be 
improved  by  decreasing  the  range  sampling  interval  but  practical  considerations  limit  the  total  physical 
number  of  range  sampling  filters  that  can  be  implemented. 
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Figure  22:  Plot  along  range  vector  containing  targets  of  AAGC  output  for  configuration  2  input  data 
with  non-deterministic  portion  of  clutter  signal  added.  Arrows  indicate  targets. 


Figure  23:  Plot  along  range  vector  containing  targets  of  DAGC  output  for  configuration  2  input  data 
with  non-deterministic  portion  of  clutter  signal  added.  Arrows  indicate  targets. 
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4.0  Summary  and  Future  Work 


In  summary,  the  use  of  the  DAGC  system  balances  the  output  response  across  the  leading  and  trailing 
edges  of  regions  with  large  cross  section  variations,  providing  a  symmetrical  response.  This  is 
accomplished  at  the  expense  of  an  increased  shadow  zone  on  the  leading  edge  of  a  sea  to  land  transition 
over  that  of  the  AAGC  response  but  with  a  significant  improvement  over  the  AAGC  trailing  edge  response. 
In  general,  the  dramatic  improvement  at  the  trailing  edge  justifies  the  small  performance  penalty  at  the 
leading  edge.  The  uniform  performance  of  the  DAGC  across  the  field  of  view  also  simplifies  interpretation 
of  the  output. 

One  of  the  major  advantages  of  a  digital  processing  system  over  the  hard  wired  AAGC  approach  is  the 
ability  of  the  DAGC  system  to  be  adjusted  to  optimize  the  output  for  any  particular  mode  of  RADAR 
surveillance.  This  ability  alone  is  strong  justification  for  considering  the  adoption  of  a  DAGC  system. 

As  discussed  in  the  report,  other  considerations  such  as  relative  stability  of  the  two  systems  and  noise 
behaviour  do  not  strongly  recommend  one  method  over  the  other.  Both  systems  should  remain  stable  over 
the  range  of  conditions  likely  to  be  encountered  during  operation  and  the  small  decrease  in  signal  to  noise 
ratio  at  the  output  of  the  DAGC  in  comparison  with  the  AAGC  system  is  likely  to  be  much  less  than  5%  for 
any  practical  implementation  utilizing  temporal  averaging  and  spatial  filtering  of  the  clutter  map. 

No  attempt  has  been  made  in  this  report  to  quantify  improvements  or  decreases  in  detection  probability  that 
would  occur  with  the  implementation  of  a  DAGC  scheme.  In  order  to  do  this  it  is  necessary  to  identify  the 
detection  scheme  being  used.  One  focus  of  future  investigations  should  be  to  quantify  the  expected 
behaviour  of  the  DAGC  system  using  different  CFAR  and  peak  detection  schemes.  In  addition,  the  use  of 
more  elaborate  statistical  noise  models  such  as  the  K-distribution  should  provide  a  clearer  picture  of 
performance  under  more  realistic  conditions.  Any  such  analysis  would  also  benefit  from  the  inclusion  of 
other  practical  receiver  dependent  parameters  such  as  bandwidth,  pulse  compression  etc. 

The  flexibility  of  the  digital  processing  scheme  also  opens  the  possibility  for  real  time  manipulation  of  the 
implemented  filter  function.  It  is  envisioned  that  the  potential  exists  to  use  edge  finding  techniques  to 
modify  the  spatial  shape  of  the  filter  interactively  to  prevent  the  filter  from  straddling  the  land-sea 
boundaries  and  thereby  further  reduce  the  shadow  zone  along  the  coast  line. 

Future  work  will  focus  on  developing  the  basic  AAGC  and  DAGC  processing  code  into  modules  which  can 
be  used  in  conjunction  with  the  DREO  GMTI  simulator  to  allow  pre-development  modeling  of  receiver 
designs  and  allow  ongoing  development  and  testing  of  detection  schemes. 
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List  of  Acronyms 


DND 

Department  of  National  Defence 

AAGC 

Analog  Automatic  Gain  Control 

AGC 

Automatic  Gain  Control 

DAGO 

Digital  Automatic  Gain  Control 

DSC 

Digital  Scan  Converter 

PRF 
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RGF 
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